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20.  Abstract  (continued) 

.j  percent  cloud  (earth)  coverage  and  cloud  heights.  Both  algorithms  were 
found  to  perform  well  under  some  meteorological  circumstances,  but  not 
so  well  under  others.  In  general,  however,  it  was  found  that  CLUSTER 
performs  better.  The  CLUSTER  algorithm  found  clear/cloud  boundaries, 
numbers  of  cloud  layers,  and  fractional  cloud  covers  that  are  closer  to 
estimates  made  by  analysts  looking  at  corresponding  IR  and  visible  images. 
The  MODAL  algorithm  found  more  cloud  layers  than  were  observed,  and 
tended  to  overanalyze  the  transition  regions  between  cloud  layers. 

Despite  the  overall  advantages  of  CLUSTER,  there  is  considerable 
room  for  improvement  and  algorithm  development  leading  to  improved 
cloud  analysis.  For  instance,  CLUSTER  finds  too  many  layeroxif  cirrus 
clouds  due  to  differing  transparencies  and  sub-pixel  size^of^irrus  cloud 
elements.  It  is  possible  to  improve  the  cirrus  Jn.yertfigproblem  by  tuning 
the  algorithm  to  be  less  sensitive  to  smaller  thermal  gradients  at  cold 
temparaTures.  CLUSTER  can  be  tuned  for  other  cloud  types  and  levels 
aawell. 

''->The  AFGWC  has  recognized  the  limitations  of  MODAL  and  the  com¬ 
plications  in  its  computer  code.  Subsequently,  a  substantially  revised 
version  of  the  3DNEPH  has  been  written,  and  is  known  as  the  Real-Time 
Nephanalysis  (RTNEPH).  With  the  continued  cooperation  of  the  AFGWC, 
we  plan  to  evaluate  this  nfew  program 's  image  analysis  algorithm  and 
recommend  changes  aniyirrynrovements  for  it,  drawing  on  our  experience 
with  MODAL  and  CLUSTERX 
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Evaluation  of  Automated  Imagery 
Analysis  Algorithms  for  Use  in  the 
Three-Dimensional  Nephanalysis 
Model  at  AFGWC 


I.  IM'RODl  C.TION 

Ail  automated  satellite  cloud  analysis  system  has  been  in  use  for  over  a 
decade  at  the  Air  Force  (ilobal  W  eather  Central  at  Offutt  \ I'  15,  Nebraska.  This 
system  is  known  as  the  M I ) N  1 : 1  ’  1 1  (Nephanalysis)  Model  and  has  been  described  bv 
Coburn  and  l  ye.-  It  is  a  very  comprehensive  system  that  ivonjes  satellite 
ini  a  eery  information  and  conventional  meteorological  data  to  produce  global 
cloud  analyses  eijlit  times  a  day.  The  infrared  processor  is  an  important  com¬ 
ponent  of  this  sy  stem,  since  the  infrared  iniaye  is  the  sole  source  of  satellite 
information  at  niuht.  Frequently  this  is  the  case  during  the  day  as  well. 

Tlie  major  component  of  the  infrared  processor  is  an  algorithm  called  MODAL 
that  was  designed  to  separate  major  regions  or  layers  in  a  2.*>  x  2a  n  mi  arid  box 
which  consists  of  an  array  of  ff  x  it  infrared  values  having  a  nominal  il  n  mi 
resolution.  The  Alol)  \  I.  algorithm  is  a  histogram  separation  procedure  that 
Uives  from  one  to  four  distinct  regions  or  layer-'  for  each  arid  box.  Our 
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objective  is  to  try  to  understand  this  algorithm  better.  What  are  its  weaknesses? 
Where  are  its  strengths? 

3  4 

\n  algorithm  developed  by  Hawkins  ’  to  analyze  satellite  data  operates  on 
infrared  images,  and  can  be  configured  to  give  much  the  same  type  of  output  as 
MODAL.  This  algorithm,  called  Cl. 1ST  III?,  separates  the  infrared  images  into 
mutually  exclusive  regions.  Ml  points  are  assigned  to  a  region,  whereas  in 
MODAL  some  points  may  be  left  unassigned.  In  CLLSTHR,  as  in  MODAL,  the 
number  of  regions  can  be  from  1  to  4.  We  will  alternately  refer  to  these  as 
layers  or  clusters.  In  MODAL  analyses  they  have  traditionally  been  called  modes 
as  well. 

The  object  of  this  report  is  to  describe  MODAL  and  CLLSTLR,  and  to  com¬ 
pare  the  results  of  .  large  sample  of  cases  where  the  two  algorithms  were  run 
side  by  side.  These  evaluations  are  made  possible  by  the  AFOL  Man-Computer 
Interactive  Data  \cces>  System  (McIDAS)  which  is  an  interactive  computer  sys¬ 
tem.  The  algorithms  were  programmed  for  the  McIDAS,  and  a  procedure  was 
arrived  at  for  comparing  the  two  algorithms  for  different  cloud  situations.  After 
some  experimentation,  a  sample  set  of  350  cases  was  selected  from  smoothed 
(  ■!  n  mi)  DMS1’  imagery.  Kach  case  consists  of  an  array  of  1G  x  16  infrared 
gras-shades,  which  means  that  each  case  contains  4  MODAL  analyses,  or  a  total 
of  1400  MOD  \[.  analyses  for  the  entire  study.  The  cases  cover  a  wide  range  of 
cloud  types,  ground  types,  and  combinations  of  cloud  and  ground.  Considerable 
insight  into  the  algorithms  was  gained  by  image  analysts  in  reviewing  these  sets 
of  data,  and  it  is  felt  that  the  basic  samples  provide  a  well-rounded  data  set  not 
only  for  this  study  but  fi  r  further  investigations. 

2.  im  \u  >  oi  mi  \i.<;oi<mi\i< 

The  basic  natures  of  the  two  algorithm.-  studied  here  are  very  different  since 
MODAL  is  a  histogram  evaluation  algorithm  and  CLLSTLR  is  a  dynamic  algorithm. 
Details  of  their  structure  and  operation  will  now  be  given. 


3.  Hawkins,  Rupert  s.  (L'T.o)  A  clustering  Technique  for  Satellite  Image 
\nal\sis,  i’roc.  flth  C  onf.  on  Weather  forecasting  and  Analysis, 

Amer.  Meteor.  Soc.^  1  1  a  -  1  ITT. 

1.  Hawkins,  Rupert  S,  (10K1)  Objective  Analysis  of  Satellite  Cloud  Imagery, 
l'roc.  lh‘11  International  Conscience  and  Remote  Sensing  Symposium, 
ILL L  i  ieoscicnce  and  Remote  Sensing  Society,  A'olume  1,  477-482. 
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2.1  Introduction  and  Details  of  MORAI. 


The  IR  processor  or  the  3D  ME  PH  merges  infrared  data,  the  surface  and  upper 

air  temperature  fields,  and  the  geography  and  terrain  fields  to  produce  an  eighth- 

2 

mesh  (25  n  mi)  cloud  analysis  on  the  3DNEPH  grid.  Cloud-top  heights,  temper¬ 
atures,  and  total  cloud  amounts  for  up  to  two  "most  significant"  layers  are 
determined.  ("Most  significant"  is  defined  in  the  3DNEPH  satellite  processor 
as  the  two  coldest  layers  separated  by  1500  feet  or  more.) 

The  automated  111  imagery  analysis  algorithm  MODA1,  is  the  routine  used 
bv  the  31) NK  I’ll.  MODAI.  is  capable  of  detecting  up  to  four  layers,  or  modes,  of 
clouds  (including  "clear"  if  it  is  present).  It  constructs  a  histogram  of  the  III 
erayshades  (see  Figure  1)  of  an  8  x  fl  III  image  array.  These  gray  shades  are 
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Figure  1.  MODAT,  Histogram  Analysis 

then  separated  into  modes  (there  are  3  modes  in  Figure  1)  and  a  representative 
grav -shade  is  determined.  I  sing  "COI.D"  as  the  coldest  grayshade  within  a  mode, 
"ML  \\"  as  the  average  grayshade  of  a  mode,  and  MAX'  as  the  most  frequent 
grayshade  within  a  mode.  Table  1  lists  the  seven  possible  ways  that  a  mode's 
representative  grayshade  is  determined  for  the  cloud/ no-cloud  decision,  and 
for  its  height,  i.acli  mode's  representative  height"  grayshade  value  is  first 
converted  to  an  equivalent  III  temperature,  then  corrected  for  water  vapor 


Table  1.  Options  for  Representative  Gray  shades  of  a  Mode 


Scheme 

Number 

Cloud/ 

No  Cloud 

Height 

1. 

COLD 

COLD 

2 . 

MAX 

MAX 

3. 

MEAN 

MEAN 

A . 

MAX 

COLD 

5. 

MEAN 

COLD 

6. 

MEAN 

MAX 

-i 

MAX 

MEAN 

attenuation  and  look  inch',  and  finally  converted  to  the  height  of  the  mode  using 
the  eighth-mesh  surface  and  upper  air  temperature  profile.  If  the  representative 
"cloud/ no-cloud "  temperature  of  a  mode  is  not  significantly  different  from  the 
surface  temperature,  then  that  mode  is  flagged  as  "no  cloud.  "  The  total  cloud 
amount  is  determined  by  summing  the  number  of  pixels  in  the  cloud  mode(s)  and 
dividing  by  the  total  number  of  pixels  in  the  image  array. 

The  choice  of  schemes  in  Table  1  allows  for  tuning  of  the  IDXfl’H  output  to 
the  needs  of  its  users.  Currently,  the  nephanalvsis  uses  Scheme  1,  Such  a 
choice  tends  to  overestimate  the  presence  of  cloud.-,  since  the  cold  point  of  a 
mode  may  be  significantly  colder  than  the  MTW  or  the  MAX  of  the  mode.  This 
scheme  also  tends  to  overestimate  the  altitude  of  low  clouds;  however,  this 
scheme  may  he  best  for  semitransparent  cirrus  clouds  since  all  of  their  gray- 
shades  may  lie  warmer  than  the  true  cloud  temperature.  An  example  is  given 
later.  \  description  of  the  logic  MODU.  uses  to  distinguish  cloud  layers  follows 
M(>I)\I.  operates  on  the  f, -hit  Ill  image  values  of  an  !!  x  H-pixel  array 
(dfi  n  mi  x  2a  n  mi  nominal  size).  Its  objective  is  to  separate  clear  regions  and 
cloud  layers  present  in  the  ima  :e  sample.  The  first  function  of  the  algorithm  Is 
to  generate  a  histogram  of  the  «  >.  :vwa\  ( see  figure  1).  Mode.-  ire  then 
defined  in  the  following  manner.  \11  searching  is  done  from  the  cold  end  of  the 
histogram  to  the  warm  end,  unless  explicitly  stated  otherwise. 

1.  \  search  i.-  made  through  the  histogram  for 

that  gray  shade  that  appears  most  frequently  in 
the  array.  This  gray  shade  is  the  mode  of  the 
mode,  henceforth  called  MW  (see  figure  1). 

The  number  of  gray  shades  of  value  MAX  must 


be  i3.  Let  Nj  =  #  of  grayshades  of  value  i. 

Thus,  if  <^'  £°  to  Step  5* 

2.  Grayshades  that  are  one  or  two  counts 
warmer  or  colder  from  MAX  are  then  un¬ 
conditionally  assigned  to  the  mode  that  MAX 
lias  begun  to  define;  however,  if  any  Niy|  \x~  l 
or  MAX ±2  =  0  (a  breakpoint),  all  further 
assignment  of  pixels  to  the  mode  on  that  side 
(warmer  or  colder)  of  MAX  is  terminated. 

3.  Further  inclusion  of  any  more  grayshades  i 
into  this  mode  farther  away  than  2  graysbade 
values  from  MAX  is  strictly  dependent  on  two 
conditions:  (1)  the  number  of  pixels  Nj  must  be 
greater  than  0,  and  (2)  the  \j  must  monotonically 
decrease  with  distance  from  M  \\.  As  soon  as 
either  of  these  two  conditions  is  not  met,  assign¬ 
ment  of  pixels  to  the  mode  of  MAX  is  terminated. 
This  is  done  on  the  warm  and  cold  sides  of  MAX 
until  either  of  the  above  conditions  is  broken  or 

t lie  end  of  the  grayshade  scale  is  reached,  which¬ 
ever  comes  first. 

4.  Go  to  Step  1,  ignoring  any  pixels  previously 
accounted  for.  1  p  to  4  modes  are  allowed.  If 

1  have  been  chosen,  go  to  Step  4. 

a.  Go  back  uui  include  any  leftover  pixels 
into  existing  mo  les  if  possible.  Leftover  pixels 
are  those  not  accounted  for  in  Steps  1-4.  Left¬ 
overs  are  placed  into  the  warmest  mode  that  is 
closest  to  it  (see  Figure  1),  and  in  no  circum¬ 
stance  will  leftover  pixels  warmer  than  the  pixels 
in  any  mode  be  included  in  a  mode.  This  ensures 
that  clear  pixels  will  not  be  inserted  into  cloudy 
modes;  however,  this  also  allows  for  some 
unassigned '  pixels. 


There  are  some  undesirable  loopholes  in  this  logic  that  are  not  immediately 
obvious.  One  situation  exists  when  several  pockets  of  IK  grayshades  with  "peaks" 
''lie  i..  lated  from  a  mode  .,o  as  not  to  be  included  in  or  assigned  to  that 
mode,  in  many  eases,  single  breakpoints  separate  substantially  large  groups  of 
pixels  from  what  should  perhaps  he  their  own  mode  (see  Figure  1)  .  The  snXLl’M 
accounts  for  these  pixels  in  total  cloud  calculations;  however,  subsequent  layer 
calculations  using  the  leftover  pixels  can  unjustly  increase  cloud  amounts  and 
heights.  Such  is  the  case  in  Figure  1,  Mode  1.  The  24 -grayshade  spread  of  this 
mode,  representing  a  T  of  approximately  40-  I  .,  is  much  too  great,  \lthough 
the  majority  of  the  pixels  in  Mode  1  are  concentrated  around  grayshade  value  .42, 
the  31 )  N  F I  ’ll  would  choose  as  a  representative  value  (using  Scheme  1  of  Table  1) 
the  coldest  grayshade  in  Mode  1,  which  is  44.  The  use  of  such  a  cold  value  for 


this  cloud  layer  seemingly  must  overestimate  its  height.  It  can  potentially  over¬ 
estimate  cloud  presence  as  well  for  modes  that  have  a  mixture  of  clear  and  cloudy 
pixels. 

I  nas  signed  pixels  appeared  in  26.7  percent  of  the  1400  MODAL  cases  we 
studied.  Most  of  the  time  the  2.T  spread  was  considerably  less  than  the  40° K  in 
Figure  1  and  many  unassigned  pixels  were  for  cirriform  clouds,  for  which  a  siz¬ 
able  ..  T  is  expected  due  to  differing  degrees  of  transparency  of  the  cirrus  clouds. 

Another  situation  that  commonly  arises  is  that  of  overanalysis.  On  account 
of  the  monotonically  decreasing  criterion  in  Step  3,  and  the  breakpoint  criterion 
in  Steps  2  and  3,  MODAL  tends  to  be  too  sensitive  to  small  features  in  the  image 
data.  Figure  1  also  demonstrates  this  point  well.  More  logically,  Modes  2  and  3 
should  be  considered  one  mode.  In  general,  MODAL  identifies  too  many  cloud 
modes  and  misses  some  clear  areas.  Some  of  this  may  be  due,  however,  to  the 
fact  that  satellite  Ill  data  is  relatively  more  sensitive  to  high  clouds  than  to  low 
clouds.  We  have  mentioned  briefly  a  few  of  the  apparent  shortcomings  of  the 
MOD  \l.  analyses.  We  realize  nonetheless  that  A  FOWL  is  aware  of  these  short¬ 
comings  and  we  anticipate  that  an  updated  III  image  analysis  algorithm  in  the 
II  I  NF  I’ll  will  attempt  to  solve  them.  (The  Ileal- 4'ime  NLl’N  IHTNFI’IL  is  a 
substantial  revision  of  the  IDNFI'll  scheduled  for  use  in  1983.) 

2.2  Tlic  Cl. I  Sill!  Msorilhm 

I’Ll  STFI1  was  designed  to  analyze  a  16  x  16  array  of  infrared  imagery.  Thus 
four  '■>  x  .'I  MODAL-type  arrays  are  contained  in  one  CI.ISTKR  array  and,  in  all 
the  algorithm  testing  we  have  done,  four  MODAI.  arrays  have  been  calculated  for 
each  t  l  .l  STI. I?  array  .  The  fact  that  ( H ‘STICK  runs  on  a  16  x  16  array  does  not 
restrict  its  potential  application  to  the  3DNFFII  since  the  critical  grayshades 
determined  by  CLFSTFH  can  be  easily  applied  to  the  four  6x8  arrays  within  the 
1  fi  x  16  array  to  give  separate  estimates  of  cloud  amounts  and  altitudes  on  the 
eighth-mesh  or  2i  n  mi  grid  scale. 

The  first  step  of  this  algorithm  is  to  separate  the  array  into  four  quadrants. 
Histograms  are  obtained  for  sixteen  intervals  per  quadrant.  The  action  of  the 
algorithm  is  to  transform  these  histograms  into  clusters.  At  each  step,  transfers 
are  made  from  quadrants  having  smallei  f mquencies  in  an  interval  to  the  quadrant 
having  the  largest  frequency  in  that  interval.  This  separates  the  data  into  unique 
clusters  completely  consistent  with  the  overall  frequency  distribution  of  the  whole 
area.  Che  critical  gray  duides,  which  are  the  transition  grayshades  between 
elll.-te  rn,  are  called  cut  levels. 

I  he  algorithm  is  essentially  sequential,  and  can  be  stated,  therefore,  as  a 
list  of  operations.  Small  loops  sometimes  occur  within  these  steps,  but  no  loop 
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occurs  among  or  between  .steps,  l'his  characteristic  makes  for  very  fast  execu 
lion  as  compared  to  algorithms  that  require  many  iterated  calculations  within 
large  loops.  The  major  steps  in  the  clustering  procedure  follow: 

STIiP  1.  Subscript  the  lfi  x  IP  infrared  array  by 
quadrants.  Phis  dimensioned  variable  could  be 
labeled  111  (  1,  (PI).  The  four  is  for  the  four  quad¬ 
rants  and  tlie  sixty-four  for  the  ft  t!  values  within 
each  quadrant. 

SPP:r  The  :1DN I. I’ll  grayshades  (fi -bit  or  HI 
Cray  shades  linear  with  temperature  from  H10  to 
din  K)  are  placed  in  Hi  intervals.  The  widths 
of  the  intervals  are  shown  in  'Pablo  1,  and  are 
recognized  to  include  all  or  nearly  all  meteor¬ 
ologically  significant  data.  Interval  1  includes 
the  warmest  counts,  and  Interval  lfi  includes 
the  coldest  counts. 

SPlil1  Pur  each  interval,  place  tin-  sum  of 
ill  four  quadrant  frequencies  into  the  quadrant 
having  the  largest  frequency  value  for  that  in¬ 
terval.  If  equals  exist  for  largest,  leave 
for  the  next  - 1 < ■  1 1.  (Vce  that  the  quadrant  data 
ire  now  being  t r.msfo rrie. !  'c  rlu-ter  data.l 

S  111’  I'.  >r  those  instances  wit  h  t  wo  or  mo  re 

1  ir.:e  a  frequencies,  place  the  stun  at  lit  it  m- 
len  il  into  the  quadrant  having  the  largest  sum 
in  udineent  intervals. 

sPIII'  a.  if  any  quadrant  (cluster)  has  less  than 
thirteen  points  (;>  percent  of  total  coverage)  place 
it  m  the  ilea  rest  cluste  r. 

si' III'  e.  I'or  any  cluster  consisting  of  two  parts 
(one  or  more  zero  intervals  between  non-zero 
frequencies),  separate  by  putting  one  in  a  zero 
cluster  if  a  zero  cluster  exists.  If  not,  put  the 
one  closest  to  another  cluster  into  that  cluster. 

NT  1.1’  7.  Combine  close  clusters.  Phis  closeness 
value  is  one  of  the  parameters  open  in  the  analysis. 

It  was  selected  so  that  clusters  consisting  of  values 
at  onl\  one  interval  and  having  a  similar  cluster  one 
interval  away  in  another  cluster  are  combined  into 
one  cluster. 

STI.’F’  ft.  eliminate  clusters  with  less  than  Hi  points 
(  i  |>ercent  of  coverage).  Phis  is  used  a  second  time 
since  the  separation  of  multiple  clusters  could  create 
a  small  cluster. 

Sl'lll’  0.  Summarize  cluster  information;  this  in¬ 
cludes  the  generation  of  histograms  of  clusters, 
number  of  clusters,  and  critical  grayshades  or 
cut  levels  that  separate  cloud  regions  from  one 
another. 

1 1 


Table  2.  Corresponding  Interval  Ranges  for  the  6-Bit  DMSP  Data  Used 
in  This  Study.  The  temperature  ranges  are  determined  by  linearly 
converting  the  ends  of  the  corresponding  grayshade  intervals  to  *K 


iterval 

DMSP  6-Bit 

IR  Grayshade  Range 

Temperature  Ran«e  °K 

16 

52-63 

227.3-210.0 

15 

49-51 

232.2-229.0 

14 

47-48 

235.4-233.8 

13 

44-46 

240. 1-237.0 

12 

42-43 

243.3-241.7 

1 1 

39-41 

248. 1-244.9 

10 

37-38 

251.2-249.7 

9 

34-36 

256.0-252.8 

8 

32-33 

259.2-257.6 

7 

29-31 

264.0-260.8 

6 

2  7-28 

267.1-265.6 

5 

24-26 

271.9-268.7 

4 

22-23 

275. 1-273.5 

3 

19-21 

279.8-276.7 

1  7-18 

283.0-281.4 

1 

0-16 

310.0-284.6 

These  cut  levels  are  the  major  result  of  the  clustering  routine.  Along  with 
the  infrared  values,  they  specify  clustered  regions  that  become  the  subject  of  in¬ 
terest,  and  represent  a  dimensional  simplification  of  the  imagery  data.  Figure  2 
shows  some  result-,  of  calculations  for  a  GOKS  III  image  sample.  Quadrant  histo¬ 
grams  appear  on  the  left,  and  cluster  histograms  on  the  right.  Summaries  are 
given  below  these  as  marked. 


SIMM  Mills  Of  I  II!  DAI  \  SKI 

The  data  set  used  for  this  study  was  taken  from  nominal  3  n  mi  resolution 
Defense  Meteorological  Satellite  Urogram  Operational  l.inescan  System  (DMSP  OI.S) 
visible  and  III  data,  only  111  data  were  actually  processed,  but  corresponding  visi¬ 
ble  images  were  used  by  the  image  analysts  to  help  in  making  accurate  subjective 
cloud  amount  and  type  classifications.  The  OI.S  is  the  primary  sensor  on  DMSP 
spacecraft.  It  is  a  dual-channel  scanning  radiometer  that  senses  reflected  light 
and  emitted  infrared  energy  in  the  0.4  to  1.0  am  and  10.2  to  12.8  am  spectral 
bands,  respectively.  A  more  advanced  design  of  the  OI.S  enables  the  earth  sampl¬ 
ing  to  vary  much  less  in  resolution  as  scanning  proceeds  from  nadir  to  higher  scan 
angles.  The  3  n  mi  oi,S  data  described  in  this  report  were  smoothed  twice:  first. 
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RANGE  «  QUAD.  FREQS.  CLSTR.  FREQS. 
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MKAN  IK 
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98.91 128.6| 179.91 

MEAN  VIS 

140.21 
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Kittlin'  llisto.jram  Data  llefore  and  \fter  Clustcrinij  With  a  Summary  lielow  of 
the  Knur  flusters.  This  example  is  from  an  K-bit  lit  DOKS  intake 


on  t!ii-  satellite  from  u.  a  n  mi  to  1 .  a  n  mi,  and  then  at  \K<!WC  from  1,5  n  rni  to 
"  n  mi.  lake  tin-  visible  data,  til.S  11!  data  are  digitized  into  one  of  lUR  possible 
j  ray.^liadi-  values  on  the  satellite,  and  truncated  to  one  of  l>4  possible  values  at 
«IWf.  Idle  II!  sensors  ire  desiijned  to  detect  a  maximum  equivalent  blackbodv 
temperature  of  10'  K  and  i  minimum  of  DO  h.  Warmer  scenes,  however,  tend 
to  h  ive  considerable  atmospheric  ittenuation  so  that  equivalent  blackbodv  tempera¬ 
ture  are  Ie  •  -  tii. oi  true  .retie  temperatures,  \ttenuation  correction  factors  vary 
f it"  1  >  it  'in  |  d> -wn  to  fat  -’1 11  I , . 

I  he  l.it a  set  eon  i  t  -  et  : .'at  <  1  I  s  Id  . If  sample.,  md  t hei r  cor respondin a  1  -1 00 
M  <  ’I )  \  1 .  .  a  - 1  pie.  taken  from  i  vide  v  arietv  <»f  atmospheric  conditions  ramtinu  from 
ole  ir  to  hiitbly  active  cloud  rejmn  .  'd.  Pile  an  effort  was  made  to  (jet  as  varied  a 

s  m  pie  as  possible,  t>  e  ei'iv  pfev  d"iii  cloud  t\pi-s  occur  niore  frei|uontly  within 
the  -et,  \1  -o,  -  i  l  ie  .  of  1  uvl-'.v  'ter  boundaries  under  clear  skies  were  selected. 

I  he  nn  |e  rl\  in^  een  -ide  r  it  mp  m  ;*  ik  in.  I  be  d  it  i  -elect  it'll  was  to  obtain  as  com  - 
I  ’  1  e 1  - ■  a  s  ip  pie  i.  pi i  -  .| b|e  f ,  ,|.  test  ii ip  1 1 1! i  ■  v i t 1  r  ■ ,  .  tuoh  as  M t  'll  \  I .  and  ( '  1 .1  ST  KR . 


» s  ' -  -  'miiii’iii'i^rii— "•  *  •  ' . 


The  samples  are  mostly  from  regions  of  feature  transition.  Flat  regions  were 
occasionally  selected  but,  since  these  do  not  generally  reveal  much  about  the 
capabilities  of  the  algorithms,  samples  that  include  more  than  one  level  or  feature 
(that  is,  land/ cloud)  wer<  ordinarily  chosen.  These  levels  range  in  contrast  from 
multilayered  cloud  scenes  to  simpler  clear/ cloud  scenes. 

The  data  set  was  sent  on  tape  by  U'tlWC  to  AFC.I.,  and  eighth-mesh  samples 
were  collected  and  saved  using  the  McIHAS.  Data  and  algorithm  results  wort- 
printed  out  for  each  of  the  3  50  eases  in  the  set.  The  results  for  Case  .102  are 
shown  in  Figure  3.  Data  such  as  shown  in  Figure  3  was  used  to  compare  tile  per¬ 
formance  of  the  algorithms.  \t  the  top  of  Figure  3  some  information  concerning 
which  imagery  the  case  came  from  is  presented.  The  data,  cloud  types  observed, 
and  background  (which  were  entered  by  the  operator  luring  the  sample  save  gen¬ 
eration)  are  given.  The  top  array  contains  the  original  fi-h'.t  infrared  image  values. 
(A  value  of  0  is  warmest,  and  (id  is  coldest.)  Ile.-ults  of  \H  >D  \  I  and  Cl. I  ST  id  H  are 
then  listed  in  coded  form.  '  I."  represents  the  lowest  layer  (warmest  mode  or 
cluster),  "  .  the  next  highest  (next  coldest)  layer,  -  the  next,  and  II  the  highest 
layer.  (Note  that  if  only  two  layer  are  found,  .  represents  the  highest  .coldest, 
one.)  i  iray  shade  ranges  for  the  four  \KMJ\I.  arrays  and  one  Cl  l  ST  1 11?  array  are 
then  given  at  the  bottom,  of  tbe  printout. 

f  igure  -i  also  illustrates  flow  the  (1.1  ST!  U  ilgonthm  would  easily  generate  a 
cloud  analysis  for  it  x  it  lit  arrays  (ei  ghth  medil  or  even  !  x  '•  arrays  (sixteenth- 
na*sh).  I  lie  critical  gray-ha  |es  or  cut  values  can  be  determined  for  the  Hi  \  Hi 
array,  as  Cl.VSn.lt  already  lee..  These  cut  values  call  then  he  applied  to  com¬ 
pute  cloud  pa  r  atr  el  e  rs  -  ueh  i  -  el  a.  i  em  ■  r  t :  ’ 1  a  •  •  -  ai  r  \  ('■  or  i  \  !  pixel 
ll beets  of  the  original  i  (>  X  Hi  nrv  r  . 

The  example  in  Figure  .3  is  for  S  made  D i’  which,  according  to  the  operator's 
notes  recorded  at  the  time  of  the  save,  consi  ;t-  of  about  a;'  percent  nirnilu  -  build¬ 
ups  and  ai)  percent  stratus  buildup.-.  I?  o  member  that  M<  >D  M.'.s  maly  -<  ■  ire  for 
four  independently  analvaed  (i  x  "  irray.-,  all  of  which  make  up  the  V>  percent  (  hi 
and  at)  percent  St,  and  are  merely  positioned  b'gethv  ”  in  Figure  :  for  co--p:iri  son 
with  the  CM  Sd'l  I!  .malvsis.  T'lis  cause  .-otre  eventual  irregularities  in  the 
layers'  representative  IK  .  ray  shade.-  !«•  occur  fro;  ;  one  \  h  hex  to  .mother, 
where. o  a  smooth  analy  sis  m  .ci;  b\  i  I  l  s  |  lit.  \ny  aio1;  -r .  dl---e  ilc 

irregularities  will  be  reflected  in  the  cloud  dialysis  output. 

d.  I  ( lllll-gorir*  el  Clear  and  Cloiiih  lb-glen. 

The  ;'.:>0  cases  were  separated  into  .-.oven  basic  categories.  Table  ,  -bows  the 
names  of  these  categories  along  with  tile  number  of  ea  .es.  From  left  to  right 
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these  are  Clear  I.and/W  ater/Snow,  Cumulus,  Stratus,  Stratocumulus,  Cirrostratus, 
Cirrus,  anti  Cumulonimbus. 


Table  a.  Composition  of  the  (  loud  Imagery  Sample  Set 
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This  separation  is  done  on  the  basis  of  the  major  classification  (the  most 
dominant  type)  for  each  sample.  Minor  constituents  usually  occur  in  all  cases 
except  clear  land  or  clear  water.  We  feel  that  this  sample  is  more  than  adequate 
for  evaluating  M ( >1  >  V I .  and  CI.l'STI  IK-type  algorithms.  The  fundamental  drawback 
in  the  evaluation  of  such  algorithms  is  the  subjective  approach  that  is  required  to 
settle  on  (a)  what  is  desired  of  an  algorithm  like  MODA1.,  and  (b)  once  decided 
upon,  how  well  MODAl.  and  other  algorithms  meet  those  criteria.  The  criteria 
we  feel  to  be  of  importance  are  enumerated  in  the  following  paragraphs.  Also 
stated  is  the  rationale  for  separating  a  grid  box  into  regions  or  layers. 

Separation  of  basic  features  for  making  simple  accurate  statements  about 
grid  boxes  is  the  goal.  In  the  iDM.l’ll  these  features  (regions)  are  referred  to 
as  layers.  Naturally,  layers  standing  out  from  other  layers  are  most  desired  to 
he  separated.  Two  layers  hardly  distinguishable  from  each  other  should  be 
analyzed  as  one  layer.  Doth  the  size  of  the  area  and  how  much  the  area  stands 
(Hit  from  adjacent  areas  should  determine  whether  or  not  it  qualifies  as  an  in¬ 
dependent  layer.  These  goals  are  not  stated  explicitly;  however,  they  are 
implicitly  included  in  algorithms  that  perform  the  separation. 

The  separation  should  he  mutually  exclusive.  In  other  words,  the  sum  of  all 
the  layers/clusters  should  constitute  the  whole  .grid  box  area.  Among  other  things, 
this  makes  for  a  complete  analysis  and  more  meaningful  data  reduction  than  if 
this  is  not  a  requirement.  PI.T'STKI?  has  this  feature.  MODAl,  lacks  this  feature 
since  sometimes  Ill  values  are  not  assigned  to  modes;  however,  subsequent  proc¬ 
essing  in  the  .‘IDN’Kl’H  decides  whether  the  unassigned  III  values  are  clear  or 
cloudy. 

The  sensitivity  of  the  separating  algorithm  should  be  reasonable  from  the 
warmest  equatorial  regions  to  the  coldest  polar  regions.  It  would  also  be  an 


lfi 


advantage  if  it  could  be  tuned  to  give  a  desired  response.  Both  MODAL  and 
CLl  ST KR  are  open  to  tuning,  but  CLUSTER  is  far  more  flexible  in  this  respect. 

For  operation  in  the  3DNEPII  model,  fast  computer  times  are  required  to 
perform  the  analysis  due  to  real-time  processing  constraints.  If  an  algorithm  is 
not  fast  enough  it  cannot  qualify  for  use  by  the  3DNKPH.  Of  course,  MODAL 
satisfies  these  requirements.  Our  experiments  with  run  times  show  that,  on  an 
average,  MODAL  takes  0.02G  seconds,  and  CLUSTER  takes  0.094  seconds.  Com¬ 
paring  area-weighted  times,  CLUSTER  runs  a  little  faster  than  MODAL  does. 

•  i.2  Iw.ilu.itiiin  of  (In-  Thu  Algoritliiu*  in  Terms  of  tlx*  Seven  H.1'11  Categories 

Clear.  From  the  point  of  view  of  region  separation,  this  is  a  simple  category. 
Problems  arise  not  from  algorithm  inadequacies,  but  from  the  very  nature  of  the 
IR  data.  There  are  instances  when  low  clouds  blend  into  the  surface  temperatures 
and,  consequently,  are  not  separated  by  either  MODAL  or  CLUSTER.  Therefore, 
clear  regions  and  scattered  cumulus  conditions  are  often  difficult  to  separate 
satisfactorily.  Figure  4  illustrates  how  the  1G  x  10  CLUSTER  array  size  has  a 
"context'  advantage  over  MODAL,  when  scattered  cumulus  are  present  in  an 
otherwise  clear  area.  CLUSTER  identifies  two  distinct  layers,  compared  to  just 
one  layer  for  MODAL. 

In  warmer  regions  where  land  temperatures  approach  those  of  water,  subtle 
land/water  distinctions  art*  barely  discernible  by  the  IR  sensor,  and  commonly 
indiscernible  by  the  MOD  A  (.and  CI.ISTER  algorithms  (Figures  and  G).  Wherever 
temperature  contrasts  are  great  enough,  neither  MODAL  nor  CLUSTER  has  serious 
problems  in  land/ water  separations. 

(  umulus.  There  are  several  problems  with  the  ability  of  both  MODAL  and 
CLUSTER  in  the  analysis  of  Cu.  Most  of  these  problems,  however,  are  due  to 
the  physical  nature  of  the  clouds  themselves.  They  generally  do  not  fill  the  field 
of  view  of  the  sensor.  Thus,  calculations  of  total  cloud  are  difficult  to  perform, 
especially  on  a  pixel -by -pixel  basis  as  the  3DNEPH  essentially  operates.  An  IR 
pixel  that  represents  a  Cu  cloud  that  partially  fills  that  pixel  contaminates  the 
grayshade:  the  temperature  that  grayshade  represents  is  actually  representative 
of  some  cloud  and  some  ground  radiation.  The  observed  temperature  is  warmer 
than  the  actual  cloud  element's  temperature,  but  colder  than  the  underlying  sur¬ 
face  temperature.  When  the  algorithms  compare  such  a  temperature  to  an  adja¬ 
cent  clear  pixel,  the  difference  in  temperature  (grayshade)  may  or  may  not  be 
enough  to  allow  MODAL  or  CLUSTER  to  distinguish  between  a  partially  cloud- 
filled  region  and  a  clear  region.  It  should  be  noted  that  this  is  not  necessarily 
an  algorithm  problem,  hut  a  physical  one  due  to  IR  sensor  characteristics,  field 
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of  view,  and  atmospheric  attenuation.  Nonetheless,  CLUSTER  seems  to  handle 
these  situations  more  favorably  than  MODAL,  since  MODAL  "tries  too  hard"  to 
separate  subtle  changes  in  grayshade.  Separation  of  Cu  from  higher  layers  does 
not  present  a  problem  for  either  algorithm. 

Figures  7  and  fi  illustrate  the  action  of  the  two  algorithms  on  cumulus  fields. 
First,  notice  that  MODAL  overanalyzes  (provides  seemingly  irrelevant  detail) 
in  some  places,  and  that  CLUSTER  gives  an  apparently  neater  analysis.  Points 
not  analyzed  into  a  level  by  MODAL  are  blank  in  the  illustrations.  These  points 
are  warmer  than  the  warmest  mode  (see  Step  a  of  MODAL  in  Section  2.  1).  While 
the  two  cases  represent  specific  situations  that  are  not  that  frequent  in  the  data 
samples,  they  do  indicate  MODAL's  apparent  problem  of  overanalyzing  the  image 
data.  MODAL  might  be  characterized  as  being  "oversensitive"  to  small  changes 
over  small  distances.  Some  of  the  variability  in  the  MODAL  analysis  occurs  at 
the  boundaries  of  the  quadrants.  MODAL  is  not  at  fault  here  since  the  individual 
quadrants  may  have  appropriate  layers.  It  must  be  remembered  that  the  MODAL 
mapping  is  made  up  of  four  independent  runs.  Symbols  may  not  represent  the  same 
temperature  range  from  one  quadrant  to  the  next. 

Stratus.  Neither  MODAL  nor  CLUSTER  has  a  problem  with  flat  fields  of  data. 
Stratus  is  one  cloud  type  on  which  the  two  algorithms  behave  similarly-.  There 
are  cases  to  lie  found  throughout  the  stratus  samples  that  support  this  observation; 
however,  again  there  is  a  tendency  for  MODAI.  to  be  oversensitive,  and  put  in  more 
layers  than  necessary. 

■Stratocumulus.  Stratocumulus  poses  less  of  a  problem  than  cumulus  since  these 
clouds  tend  to  fill  more  of  the  field  of  view.  Thus,  many  of  the  physical  con¬ 
straints  inherent  in  cumulus  detection  and  analysis  are  not  so  prevalent  for  the 
proper  detection  and  analysis  of  Sc.  Doth  MODAL  and  CLUSTER  tend  to  find  too 
many  layers  of  Sc.  They  do  not  always  do  so  on  the  same  cases.  Examples  of 
MODAL  over.uialy sis  of  Sc  are  in  Figures  9,  10,  and  11,  and  for  CLUSTER  in 
Figures  12,  12,  and  1-1.  Doth  are  finding  subtle  changes  in  the  temperatures  of 
cloud  fields.  We  feel  that  the  cloud  analysis  would  be  improved  if  this  sensitivity 
was  tuned  out  of  the  algorithm  to  be  used  in  the  3DNICPH. 

Oir rostratus.  Cirro-stratus  fields  tend  to  appear  flat  and  are  generally  well- 
analyzed,  since  temperature  differences  among  these  ice  clouds  and  most  any 
other  categories  in  an  image  are  great.  While  both  algorithms  perform  well, 
MODAL  has  the  tendency  to  find  too  many  layers,  as  noted  in  other  categories 
(see  Figure  IS,  for  example). 
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MOD\l.  and  CU  STKK  Analysis  for  Case  199 


Cirrus.  Both  algorithms  tend  to  find  too  many  layers  of  cirrus.  MODAI.'s  over- 
analysis  appears  to  be  worse  than  CLlSTER's,  however.  Cirrus  boundaries 
naturally  are  the  cause  for  this  problem,  as  transition  zones  among  solid  cirrus 
and  other  regions  are  partially  filling  the  sensor  field-of-views.  Also,  cirrus 
variations  in  thickness  and,  hence,  radiative  temperature  contribute  to  this 
problem.  While  some  altitude  variation  of  cirrus  is  expected,  the  algorithms 
respond  too  strongly  to  field-of-view  and  emissivity  variations.  As  with  previous 
cases,  these  sensitivities  should  be  tuned  out  of  the  operational  algorithm. 

Cumulonimbus.  Both  algorithms  have  little  problem  detecting  the  presence  of 
cold  regions  adjacent  to  warm  regions;  however,  both  algorithms  tend  to  be  too 
sensitive  to  Cb  edges.  MODAL  seems  to  find  more  layers  at  the  C'b  edges 
(Figure  Hi)  than  CI.l  STLR.  It  is  interesting  to  note  that  in  many  Cb  cases  where 
MOD  \L  has  "undefined"  pixels  (see  Step  f>  of  MODAL  in  Section  2,  1),  clouds 
are  generally  found.  In  addition,  these  areas  seem  to  be  where  CI.ISTLR 
ove  ranaly  zes.  ’ 

■  Li  t‘'m|iH*iirirs  til  ;in<l  'lodes 

The  i  >0  cases  were  reexamined  to  study  the  tendency  of  MODAL  and  CI.ISTLR 
to  find  too  many  cloud  layers.  For  this  study,  we  looked  at  clear  cases  that  were 
entirely  clear,  anti  cloudy  cases  with  only  one  category  of  cloud.  We  were  left 
with  4  7  clear  cases;  27  with  a  mixture  of  land  and  water,  and  20  with  a  uniform 
background.  We  also  found  l(iH  suitable  cloudy  cases.  The  cloudy  cases  could  be 
overcast  or  partly  cloudy,  as  long  as  the  background  was  uniform  and  only  one 
cloud  category  was  present  in  each  case.  The  objective  was  to  identify  cases 
where  the  MclD  \S  observer,  using  both  IK  and  visible  imagery,  found  only  one  or 
two  regions  ptr  case,  so  that  only  one  or  two  clusters  or  modes  would  be  the 
desired  result  of  the  algorithms. 

Figure  17  shows  the  percent  distributions  of  the  numbers  of  clusters  and 
modes  for  these  cases.  Cirro-stratus  and  cumulonimbus  are  not  shown  since 
there  were  too  few  cases  to  define  a  distribution.  CI.ISTLR  performed  well  for 
cumulus  cases  since  the  desired  one  or  two  clusters  were  found  8G  percent  of 
the  time.  For  clear  cases,  both  MODAL  and  CLFSTFK  performed  well.  There 
were  a  few  cases  over  cold  land  or  snow  where  CLFSTFR  found  two  or  three 
clusters  when  only  one  was  desired.  The  larger  array  size  for  CI.I  STLR  may  be 
more  sensitive  to  thermal  gradients  over  clear  land  than  the  smaller  array  for 
Mol)  \L.  We  will  look  at  this  in  future  studies  to  ensure  that  the  colder  land 
clusters  are  not  improperly  identified  as  clouds. 
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fijMBER  OF  CLUSTERS.' MODES  CLUSTER  MODAL 


I-' i  jure  17.  Percent  I  »istribution  of  the  Numbers  of  (  lusters/  Modes  for  Entirely 
Clear  I  ast’s  and  Cloudy  Cases  W  ith  <  >nlv  ( >ne  (  aterorv  of  (  loud  Per  Case, 


The  cloudy  categories  in  Figure  1?  confirm  the  examples  in  the  previous 
section.  Moth  MODAL  and  (’Ll 'STICK  are  tending  to  overanalyze  cloud  layers, 
finding  11  or  -4  modes  or  clusters  when  only  1  or  2  are  desired.  Moreover,  MODAL 
is  overanalyzing  more  frequently  than  OLl'STKR.  The  cirrus  cases  are  the  worst 
(The  few  cases  of  pure  cirrostratus  and  cumulonimbus  also  tended  to  have  2  or  4 
modes  or  clusters.)  It  can  be  argued  that  cirrus  may  be  found  in  multiple  layers 
so  that  3  or  4  clusters  or  modes  would  be  acceptable  for  some  cases.  Nothing  in 
our  meteorological  experience,  however,  leads  us  to  believe  that  this  should 
happen  !'2  percent  of  the  time  as  Figure  17  shows  when  the  percentages  for  3  anti 
4  clusters/ modes  are  summed. 

Two  further  comments  related  to  the  excessive  number  of  modes  found  bv 
MOD  \[,  should  lie  made.  First,  the  most  recent  technical  report  on  the  :{DNI!1’II“ 
states  that  a  mode  is  defined  for  a  single  grayshade  that  contains  n  or  more 
samples,  or  a  group  of  adjacent  gray  shades  containing  12  or  more  samples.  This 
is  no  longer  the  case  in  M('I)\L,  since  a  single  grayshade  that  contains  3  or  more 
samples  can  presently  define  a  mode.  These  "small”  modes  contribute  heavilv  to 
Ml 'DAI. 's  tendency  to  find  more  modes  than  observers  can  find  in  the  111  or  visible 
imagery.  The  second  comment  is  that  not  all  the  modes  in  MODAL  are  designated 
as  separate  cloud  layers  by  the  IDM.l’II.  There  are  somewhat  arbitrary  rules 
elsewhere  in  the  IDNLI'II  code  for  combining  the  mode.-,  into  no  more  than  two 
cloudy  layer  -.  Although  wo  did  not  duplicate  these  rules  in  ibis  study  ,  we  suspect 
that  the  overunalysis  of  modes  will  still  be  detrimental  to  the  final  cloud  layers  in 
the  3D  \ I : I’ll,  at  least  when  cirrus  in  present,  and  prohabl-  for  other  cloud  cate¬ 
gories  as  well. 

\ s  anticipated,  both  algorithms  performed  extremely  well  (over  ho  percent 
accuracy)  when  analyzing  samples  where  only  one  mode  or  cluster  was  desired. 

It  should  lie  noted,  however,  that  the  majority  of  eases  containing  only  one  homo¬ 
geneous  region  were  clear  land,  water,  or  snow  and  uniform  stratus,  imagery- 
types  tit  at  present  only  minor  problems  to  the  algorithms. 

4.1  I’ri'liinin.iri  Lloml  Lmi-r  Ix.ilmilioti 

While  cases  were  being  identified  and  saved  on  the  Moll)  \S,  the  image 
analysts  made  subjective  estimates  of  the  fraction  of  cloud  categories  over  the 
Ifi  x  If;  array  of  IK  values.  The  observer  could  look  at  both  visible  and  II!  images 
so  that  the  estimates  of  cloud  cover  for  low  clouds  were  generally  trustworthy , 
and  were  assumed  accurate  to  tJO  percent.  \s  in  tile  previous  section,  simple 
cloud  cases  with  only  one  cloud  category  appearing  in  each  ease  were  studied. 

An  observer  reviewed  each  case  for  both  M<  >D \l.  and  I'Ll  STL  II,  identified  the 
cloudy  modes  and  clusters,  and  calculated  the  total  cloud  cover  for  each  1  f>  x  lfi 


M 


array  baaed  on  the  numbers  of  samples  in  the  cloudy  modes  and  clusters.  The 
observer  was  performing  the  role  of  the  3DNHPI1  surface  temperature  field  in 
separating  the  clear  and  cloudy  modes  or  clusters.  (This  was  a  time-consuming 
technique  that  we  will  endeavor  to  improve  for  future  cloud  analysis  studies  on 
the  MclDAS. )  Results  are  given  in  Tables  4  and  a. 


Table  4.  Cloud  Cover  Within  10  Percent  of  Observed 


Table  .">.  Cloud  Cover  Closest  to  Observed 


Table  1  shows  the  percentages  of  cases  that  had  cloud  cover  within  10  percent 
of  the  observed.  For  the  total  of  the  four  categories  of  clouds,  CLUSTER  per¬ 
formed  better  than  WOD\L  (.">3  percent  compared  to  4.3  percent)  but  there  is  ob¬ 
viously  room  for  improvement  for  both  algorithms.  The  cloud  cover  estimates 
were  weakest  for  cumulus  and  cirrus  categories,  which  is  expected  since  these 
categories  are  most  likely  to  have  subpixel  cloud  elements  and  highly  variable 
changes  in  cloud  emissivities  from  one  pixel  to  the  next.  The  categories  with 


larger  cloud  elements,  St  and  Sc,  had  more  accurate  cloud  cover  estimates  for 
hot!)  algorithms. 

I'alile  .I  shows  the  percentages  of  cloud  cover  estimates  that  were  closest  to 
the  observed  values  on  a  case-by-case  basis.  Cl.l'STI.Ii  liad  cloud  cover  esti¬ 
mates  closest  to  observed  for  3  out  of  t  categories,  and  for  the  sum  of  all  cate¬ 
gories;  however,  the  differences  between  M('l).\l.  and  01.1  Si'Ll?  are  not  great 
in  this  comparison. 

The  average  cloud  covers  for  observed,  MODAL,  and  Cl.TSTLH  estimates 
were  calculated  for  each  category  to  sec  if  the  algorithms  were  finding  more  or 
loss  clouds  than  the  observed.  Very  little  difference  was  found  for  Lu,  St,  and 
Sc,  but  \loi)\l.  mid  Cl. I  Si'Ll!  both  found  about  ■<  percent  more  (  i  than  the  image 
observers,  I'he  image  observer  estimates  of  Ci  coverage  often  seemed  low  com¬ 
pared  to  wind-drawn  contour  analyses  of  the  Hi  x  Hi  arrays  of  111  data.  For  this 
cloud  category,  moreover,  the  visible  data  was  not  very  helpful  for  estimates 
of  cloud  cover  from  the  in:  ige  displays. 

file  siihiective  estimates  of  percent  cloud  cover  were  almost  indispensable 
for  evalu  iting  the  performance  of  I'Ll  Si'Ll!  and  MODAL,  and  have  convinced  us 
that  inter  ictiw  computer  procedures  should  he  established  for  improved  cloud 
cover  verification  in  future  studies. 


I  COM  I  l  'IONS 


\\  e  have  evaluated  the  iDXLIMI  MOD  H.  algorithm,  along  with  an  alternative 
LI. I  SILK  algorithm,  and  have  found  that  Cl. I  Si'Ll!  performs  better.  The 
(’Ll  Si'Ll!  algorithm,  finds  clear/ cloud  boundaries,  numbers  of  cloud  layers,  and 
fractions  of  cloud  cover  that  are  closer  to  those  estimated  b\  observers  evaluating 
the  II!  and  visible  images.  I' he  Mi  *i)  V I .  algorithm  finds  more  cloud  layers  than 
arc  observed,  uni  tends  to  overanalyze  regions  of  transition  from  one  layer  to 
another.  The  I'Ll  Si'Ll!  algorithm  is  designed  to  assign  all  points  to  a  layer  in 
an  optimal  manner.  Lite  M(>I)\L  algorithm  frequently  (27  percent  of  the  time) 
terminates  its  histogram  mahsis  with  unassigned  points.  These  points  are  later 
assigned  to  warmer  modes,  which  will  tend  to  overestimate  the  altitudes  of  clouds 
in  those  modes  (except  eirriform  c  itegories),  since  the  dDNI'.Fil  uses  the  coldest 
grayshade  in  a  mode  to  represent  the  entire  mode.  ('Ll 'ST till  uses  a  larger  array 
of  satellite  data,  anil  appears  to  benefit  from  the  extra  information  by  improved 
detection  of  low  clouds  such  as  cumulus.  CH  STL’H,  moreover,  produces  all 
the  information  necessary  for  cloud  cover  estimates  at  the  same  mesh  size  as 
the  sDXLl’ll.  On  the  McID  \S,  CLLSTFH  takes  very  nearly  the  same  computer 
time  as  MOD  \L,  so  that  its  computational  requirements  are  reasonable. 


:!(; 


Despite  the  overall  advantages  of  C'f.l  STKH,  there  is  considerable  room  for 
improvement,  and  algorithm  development  leading  to  accurate  cloud  analysis  should 
continue.  In  particular,  CLI  STIill  finds  too  many  layers  of  cirrus  clouds,  and 
the  fractional  coverage  of  both  cumulus  and  cirrus  clouds  usually  differs  from  the 
apparent  coverage  on  the  satellite  images.  When  thermal  gradients  are  found  in 
clear  land  or  snow,  it  can  also  find  two  or  three  clusters  when  only  one  is  desired. 
The  MODAL  analyses  were  slightly  better  in  these  areas. 

I  nlike  MODAL,  the  CKFSTKH  algorithm  can  be  tuned  in  a  number  of  places 
to  reduce  sensitivity,  and  give  results  closer  to  those  that  actually  exist.  For 
instance,  it  might  lie  possible  to  improve  cirrus  layering  by  increasing  the  in¬ 
terval  sizes  in  the  cold  temperatures.  Surface  regions  can  be  simplified;  that  is, 
analyzed  into  fewer  clusters,  by  increasing  the  minimum  size  of  a  cluster.  These 
are  but  two  examples  of  how  CLI'STKH  is  open  to  tuning  by  making  very  simple 
alterations  to  the  algorithm.  There  are  other  places  in  the  algorithm  code  where 
simple  modifications  can  be  made  to  better  model  specific  meteorological  condi¬ 
tions  ttiut  occur  commonly  on  a  global  basis. 

The  M’CWC  has  recognized  the  limitations  of  MOD\L,  and  the  complexity  in 
computer  code,  so  that  a  substantially  revised  version  lias  been  written  for  the 
planned  upgrade  of  the  HDXKl’M  known  as  the  Heal -Time  Nephanalysis  (H  TMil’II). 

With  the  continued  cooperation  of  XFLWC,  and  drawing  on  our  experience  with 
MODAL  and  LLLSTL'K,  we  plan  to  evaluate  the  IM'M.l’II  satellite  analysis  algorit.hu 
We  also  plan  to  use  the  capabilities  of  the  McIDAS  system  for  improving  estimates 
of  cloud  truth  for  the  evaluation  and  refinement  of  the  cloud  analysis  algorithms. 
Finally,  we  shall  propose  the  use  of  a  modified  ('Ll  STKH  algorithm,  developed  at 
\F(iF,  in  replacement  of,  or  in  tandem  with,  the  H  l'\ il  1 M 1/ HDN Fi Ml  image  analysis 
algorithms. 
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